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THE  ATOMIZATION  OF  LIQUID 
AND  TWO-PHASE  FLOWS1 


M.  S.  Volynskiy 


In  chambers  with  gas  flow,  mixing  in  the  zone  of  motion  and 
evaporation  of  the  drops  precedes  the  combustion  of  the  liquid  fuel 

In  this  article  we  mainly  examine  mixing  under  conditions  of 
supersonic  flows.  A  number  of  the  devices  (certain  chambers,  a 
jet  engine  nozzle)  is  characterized  by  a  two-phase  mixture  flow 
under  operating  conditions  at  a  Mach  number  of  M  >  1  [l].  This 
phenomenon  which  consists  of  the  atomization,  the  motion  and  the 
vaporization  of  drops  in  a  fuel  flame  and  in  which  thermal,  gas- 
dynamic  and  aeromechanical  factors  are  combined,  requires  a  complex 
investigational  method. 


In  this  article  the  author  mainly  examines  the  physical  side 
of  the  problem  -  the  construction  of  a  model  of  the  phenomenon, 
making  it  possible  to  reduce  the  real  flow  to  an  idealized  two- 
phase  gas-drop  flow  setup.  The  process  being  investigated  is 
broken  down  into  its  basic  stages  described  by  its  own  systems  of 
equations  which  are  solved  on  computers.  The  results  of  the  cal¬ 
culations  are  compared  with  the  experimental  data  on  the  injection 
of  a  liquid  into  a  free  jet  having  supersonic  velocity. 


!V.  N.  Suchkov  and  R.  I. 
part  of  the  work. 


Piskun  took  part  in  the  experimental 
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Designations 


m',  w\ 


y,  yB,  z,  zB  -  the  coordinates  of  the  flame  boundaries 
and  of  the  shock  wave; 

t  —  the  time  of  drop  motion; 

y#’  Z#  "  p?aene°sdJSfaendOx02f  flame  asymPtote 

W*  P*  P’  T*  M*  °P  '  t^°C^y*  P^ssure,  density,  tempera¬ 
ture,  Mach  number,  specific  heat  of  the 
gas; 

T"  M|’ rpl' lq  ~  velo°ity,  pressure,  density,  the  tempera- 
Mach  number,  specific  heat,  P 

index  in  the  incoming  gas  flow 
(initial  cross  section); 

Vl*  al*  V  c«.  -  velocity,  diameter,  temperature 

1  1  specific  heat  of  a  drop;  *mPerature, 

V»  a*  Tw*  V  P*  “  the  current  velocity,  diameter  and 

temperature,  the  specific  heat  and 
density  of  a  drop; 

4P  "  feed^reS8Ur6  dlfferentlal  In  the  liquid 

S*  -  the  concentration  of  liquid  in  the 
asymptotic  section  of  the  cold  flame; 

dc  -  the  diameter  of  the  atomizer  aperture; 

1  '  V'  P’  Ar‘  k'  ~  the  flow  rate,  velocity,  temperature 

Son,SfUr\?enSity*  Mach  number  of  gas, 
in  She  llquld*  a«l«l>ntlc  index 

part  of "oaf  section  of  the  main 
part  of  the  atomization  flame; 

p0»  TQ  -  pressure  and  stagnation  temperature  in 
the  incoming  flow; 

X  ~  tiorr lat6d  interval  of  total  vaporiza- 

o  -  surface  tension  of  the  liquid; 

L3  "  «SrntSl  tnterval  of  total  vapori- 

Z  -  the  volatility  of  the  liquid  (fraction 

Of  vaporized  substance);  ^ctlon 

m#  ~  sectiinWnfaJh  ln  ?Je  asyraPtotic  cross 
section  of  the  cold  flame; 

m»  -  gas  and  liquid  flow  rates; 

Mn,c  ~  shock; Umber  °f  the  SaS  bey°nd  the  normal 
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We  -  critical  Weber  number. 

«p 


INJECTION  OP  A  LIQUID  INTO 
A  SUPERSONIC  PLOW 

Work  [1]  set  forth  the  experimental  procedure  and  certain 
results  of  the  investigation  of  the  injection  of  a  liquid  into  a 
supersonic  flow  in  the  absence  of  phase  transitions  (vaporization). 
The  liquid  is  fed  through  a  direct- jet  aperture  (with  the  center 
at  the  origin  of  the  coordinates)  from  the  surface  of  a  well  stream¬ 
lined  atomizer  along  axis  Oy,  at  an  angle  90°  to  the  flow  direction 
(axis  Ox).  Here  the  zone  of  boundary  layer  separation  near  the 
root  of  the  atomization  flame  is  small  and  the  shock  wave  is  close 
to  a  normal  shock.  Photographs  (made  by  Toepler’s  method)  show  the 
presence  (Fig.  1)  of  an  external  flame  boundary  which  is  sharply 
outlined  and  approaches  the  horizontal  asymptote. 

The  processing  of  the  photographs  makes  it  possible  to  plot 
dependences  for  the  ordinates  of  the  flame  asymptote  y*,  z*  in  two 
projections 

y'ldc\  t’ldc  -  A  (ZHP/Wllr  K\  ( 1 ) 

in  the  parentheses  -  the  ratio  of  the  ram  pressure  of  the  liquid 
stream  and  of  the  gasj  A  and  n  -  experimental  constants,  different 
for  planes  xOy  and  xOz. 

Then,  the  initial  sections  of  the  flame  boundaries  with  the 
coordinates  y,  z  and  of  the  shock  wave  y_,  z_  can  be  described  by 
the  approximate  criterial  relationships 


yly,  yjv‘  «  *'(*//;  m,) 

(2) 

s/x*.  t,Jz‘  «  <P  (rlz'\  M,) 

(3) 

within  the  limits  0  <  x/y*;  x/z#  <.  2-3.  Figure  2  shows  a  curve  in 
dimensionless  coordinates  yfl/y*;  x/y*  for  M1  ■  2  and  2.9. 
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GRAPHICS  NOT  REPRODUCIBLE 


Pig.  1.  An  alcohol  atomization  flame  in  a 
supersonic  air  flow  (M  =  2.9,  TQ  =  293°K, 

PQ  =  16  atm(abs . ) ,  AP  =  40  atm(gage),  dc  = 

=  1  mm) . 


*/>' 


Fig.  2.  Shock  wave  of  the  atomization  flame  (TQ  =  293°K). 
1  -  alcohol  (M^  =  2.9,  AP  =  40  atm(gage),  dc  =  1  mm,  PQ  = 

=  16  atm(abs . ) ) ;  2  -  alcohol  (1^  =  2.9,  AP  =  15  atm(gage); 
d£  =  1  mm,  PQ  =  16  atm(abs.));  3  -  gasoline  (1^  =  2.9, 

AP  =  23  atm(gage),  dc  =  2.5  mm,  PQ  =  16  atm(abs.));  4  - 
gasoline  (M1  =  2.9,  AP  =  40  atm(gage),  dc  =  0.4  mm,  PQ  = 

=  16  atm(abs . ) ) ;  5  -  water  (M1  =  2.85,  AP  =  40  atm(gage), 
dc  =  1  P0  =  ^7  atm(abs.));  6  -  water  (K^  =  2.85,  AP  = 

=  65  atm(gage),  dR  =  1  mm,  PQ  =  47  atm(abs.));  7  -  alcohol 
(M1  =  2,  AP  =  40  atm (gage ) ,  dc  =  1  mm,  PQ  =  12  atm(abs.)); 
8  -  alcohol  (M^  =  2,  AP  =  20  atm(gage),  dQ  =  1  mm,  PQ  = 

=  12  atm(abs . ) ) . 
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Work  [2]  has  examples  of  dependences  of  the  type  In  (2).  An 
extremely  Important  characteristic  of  the  injection  process  is  the 
coarseness  of  the  atomization  (appropriate  formulas  are  given  in 
work  [1]).  The  coarseness  of  the  drops  (5-15  pm)  is  substantially 
less  than  the  coarseness  of  the  particles  in  subsonic  atomization 
(diameter  20-150  pm).1  Knowledge  of  the  flame  shape,  the  shock 
wave  and  drop  coarseness  will  make  it  possible  to  set  up  a  method 
for  rating  the  flow  within  the  flame  in  the  case  of  phase  transi¬ 
tions. 

THE  THEORY  OP  A  TWO- PHASE 
MIXTURE  PLOW 

The  process  in  a  certain  (main)  section  of  the  atomization 
flame  in  a  heated  flow  can  be  approximately  compared  to  a  steady- 
state  one-dimensional  flow,  in  which  occur  phase  transitions 
(vaporization).  Plows  of  this  type  are  encountered  in  various 
technical  devices.  The  investigation  of  a  two-phase  gas-drop  or 
gas-solid  particle  mixture  is  of  interest  per  se, 

In  contrast  to  the  original  setup,  the  simplest  setup  for 
calculating  the  vaporziation  of  a  single  drop  [2]  in  this  type  of 
setup  takes  into  account  not  only  effect  of  the  gas  on  the  drops, 
but  also  the  effect  of  the  drop  on  the  gas,  which  leads  to  a  vari¬ 
ation  in  its  parameters:  velocity,  temperature,  pressure,  the 
composition,  and  consequently  the  rate  of  evaporation.  The  depen¬ 
dence  of  the  physical  constants  on  the  composition  and  the  tempera¬ 
ture  of  the  gas  are  accordingly  considered.2 

The  drops  do  not  interact  with  each  other.  The  mixture,  for 
simplicity,  is  taken  to  be  monodisperse;  the  diameter  spectrum 
can  be  considered  in  the  general  case.  The  gas  and  drop  flow  with 


‘The  discussion  concerns  an  average,  median  diameter  of  drop 
spectrum. 

2 A  similar  approach  was  employed  by  G.  A.  Varshavskiy. 
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phase  transitions  is  described  by  a  system  of  three  differential 
equations:  equations  of  motion,  heating  and  drop  vaporization  [2]. 
and  conventional  equations  of  one-dimensional  steady-state  gas 
dynamics:  conservation  of  energy  consumption,  of  the  momentum 
of  a  gas-drop  mixture  and  gas  state  equation. 

The  parameters  of  the  drops  a,  v,  Tw  and  of  the  gas  w, ,  T  ,  p 
will  be  the  initial  and  boundary  conditions  in  an  initial  cross*  1 
section  of  area  F±  and  the  shape  of  the  channel  is  P  -  P(x).  The 
desired  values  are:  volatility  Z  -  the  portion  of  evaporated  sub¬ 
stance,  characteristics  a,  v,  T^,  w,  T,  P  as  a  function  of  coordi¬ 
nate  x  or  time  x  and  the  interval  of  total  vaporization  L  when  z  - 
=  1  and  a  *  0. 


Some  results  of  the  solution  of  this  system  for  flow  in  a 
cylindrical  pipe,  obtained  with  computers,  are  depicted  on  graphs 
(Pigs.  3  and  4).  A  two-phase  mixture  consists  of  an  air  stream  con¬ 
taining  drops  of  gasoline.  The  obtained  dependences  are  conven¬ 
iently  presented  as  functions  of  dimensionless  time 

t  =  t/t*,  t*  =  af/Oi,  (1,1) 

where  x  -  the  time  of  the  motion  of  a  drop;  x«  -  the  time  of  total 
static  vaporization  of  a  drop  with  a  diameter  of  a,  at  equilibrium 
temperature  T^  =  Tp_^  (Tp^  -  the  temperature  of  isothermal  vapori¬ 
zation  of  the  liquid  which  corresponds  to  the  initial  cross  section) 
using  invariable  Initial  gas  parameters. 

Figure  3  presents  the  volatility  curves  Z;  the  drop  diameters 
are  characteristic  for  atomization  in  a  supersonic  flow. 


1 Function  $1  is  readily  calculated  from  known  formulas  of  the 
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Pig.  3*  The  volatility  of  gaso¬ 
line  in  a  supersonic  flow  (g  = 

■  1.3,  M  *  2.8,  Tx  -  1050°K, 

w1  ■  1780  m/s,  *  1  atm(abs.), 

Tm  -  293°K). 

1  and  4  -  a  3  16  ym;  2  and  5  - 
a  ■  10  ym;  3  and  6  -  a  =  5  ym. 
Designation:  cbh  3  s. 


The  time  of  total  vaporization  is  fractions  of  a  microsecond 
and  substantially  increases  with  an  increase  in  particle  dimensions. 
Curves  are  given  for  the  two  extreme  values  of  initial  relative 
drop  velocity  3  v1  -  w^,  namely  u1  3  0  (broken  lines)  and  u1  3 
3  620  m/s  (solid  lines).  The  first  group  of  curves  corresponds  to 
the  absence  of  the  drop  blowing  in  the  initial  section,  the  second  - 
to  initial  blowing  at  transonic  speed.  The  numerical  solution  of 
the  equation  of  motion  of  a  single  drop  in  the  parameter  field 
beyond  the  shock  wave  gives  the  relative  velocities  less  than  the 
speed  of  sound.  In  a  two-phase  mixture  flow  in  a  cylindrical  pipe 
the  relative  drop  velocities  also  nowhere  exceed  sonic  speed. 

Figure  4  gives  the  calculated  dependences  for  the  relative  gas 
parameters:  velocities  w/w1i  pressures  P/P1;  temperatures  T/^  in 
the  cross  section,  where  vaporization  is  carried  out,  and  the  length 
of  the  total  interval  X,  on  the  gas  Mach  number  M1  in  the  initial 
flow  section. 
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JiJ*  J*  Qfs-dynamic  characteristics  and  in¬ 
terval  of  total  vaporization  in  a  two-phase 
mixture  (ax  -  16  jim,  B  ■  1.8,  P1  -  l  atm(abs.) 
Ti  B  1700°K,  Twl 
air. 


293°K,  =  o);  gasoline  in 


Upon  the  transition  of  value  Mj  through  unity  an  abrupt  change 
is  detected  in  the  values  and  character  of  the  regularities  for  the 
indicated  parameters.  When  <  1  the  motion  of  the  gas  and  the 
particles  Is  decelerated  (w  <  w^,  and  the  temperature  increases, 
n  ncrease  of  the  time  the  drops  remain  In  the  heated  now  gives 
small  Intervals  X.  In  this  case  the  heavier  drops  overtake  the  gas 
°r  M1  >  1  the  flo»  velocity  Is  substantially  greater,  and  the 
temperature  less  than  In  the  first  case,  which  sharply  Increases  X. 

ATOMIZATION  FLAME  IN  A  SUPERSONIC  FLOW 

Let  us  now  move  on  to  an  analysis  of  a  real  atomization  flame 
breaking  it  down  into  three  sections:  initial,  main  and  mixed 
(Pig.  5).  In  the  initial  section  (with  a  length  of  the  order  of 
ordinate  y«)  decay  of  the  liquid  stream,  and  breaking  up  of  the 
drops  occur  which  then  move  along  disordered  and  unsteady  trajec¬ 
tories.  Vaporization  here  is  insignificant,  since  the  drop  concen¬ 
tration  is  great,  and  their  temperature  still  low.  Toward  the  end 
of  the  section  (section  F»)  the  drops  acquire  a  velocity,  close  to 
the  flow  velocity,  emerging  in  a  system  of  ordered  trajectories 
The  subsequent,  main  section  is  characterized  by  a  more  uniform 
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Fig.  5.  Plow  setup  in  an  atomization 
flame.  1  -  initial  flame  section;  2  - 
main  flame  section;  3  -  mixed  section;  KB  - 
cross  section  P^;  RN  -  cross  section;  F'  - 

the  exit  cross  section  for  the  initial 
section  and  entry  cross  section  for  the  main 
section;  MN  —  the  shock  wave  cross  section 
Fm;  CD  -  asymptotic  section  F«;  KSR  -  the 

flow  boundary  I  and  II. 

field  of  distribution  of  all  parameters  and  regular  particle  motion. 
In  this  case,  intense  vaporization  occurs  accompanied  by  cooling 
of  the  gas.  Experiment  shows  that  the  vaporization  section  is  not 
very  extended. 

We  will  therefore  assume  that  in  the  vaporization  section  the 
interaction  of  the  intra-flame  two-phase  flow  with  the  surrounding 
medium  does  not  play  a  decisive  role  and  occurs  mainly  on  the 
periphery  of  the  flame.  Employing  gas-dynamic  terminology,  the 
initial  section  can  be  compared  to  a  mechanical  nozzle,  and  the 
main  section  -  to  a  combination  of  mechanical,  thermal  and  flow 
rate  nozzle.  The  third  section  (it  is  not  examined  in  the  article) 
is  characterized  by  mixing  of  the  intra-flame  flow  with  the 
surrounding  medium. 


< 
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For  calculating  the  processes  in  the  main  section  we  use  the 
one-dimensional  theory  of  two-phase  flow  mentioned  above.  However, 
it  is  necessary  to  know  the  parameters  of  the  gas  and  drop  mixture 
at  the  entry  to  the  main  section  F*.  They  should  differ  signifi¬ 
cantly  from  the  parameters  of  the  incoming  flow. 


The  problem  of  determining  values  w',  T',  P‘,  8*  at  exit  from 
the  initial  section  arises.  Let  us  visualize  atomization  flame  as 
a  certain  "liquid  lattice”  possessing  hydraulic  resistance  which 
consists  of  wave  resistance  and  resistance  connected  with  the  flow 
of  a  drop  cluster  after  the  wave.  The  gas  flowing  from  section 
AB,  is  divided  into  two  flows.  Flow  I  penetrates  inside  the  flame, 
and  flow  II  around  about  its  boundaries  from  without.  Of  the 
entire  maximum  flow  rate  F'w^  going  towards  the  flame,  only  the 
part  m'  *  F'w'p'  penetrates  inside  it,  and  that  is  why  the  drop 
concentration  in  section  F’ 


a,  _  + 

r 

ftl 

will  be  greater  than  nominal 

Q* +% 

P - m' - • 


(5) 

(6) 


and  parameters  w',  T',  P'  -  significantly  differ  from  the  parameters 
wl»  Ti»  pi* 


Let  us  make  the  following  hypotheses  about  the  initial  flame 
section. 


1.  We  disregard  the  processes  of  phase  transition  (vapori¬ 
zation)  in  small  section  X’  and  let  us  not  consider  the  interaction 
of  flows  I  and  II. 

2.  Let  us  disregard  in  section  F'  the  vertical  components  of 
rates  w'  and  v’  which  are  small  here. 

v  J 
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3.  Let  us  not  consider  the  frictional  forces  on  the  monitoring 
surfaces  and  let  us  disregard  the  reduction  in  the  cross  sectional 
flow  area  F'  caused  by  the  drops. 

The  selection  of  length  X'  is  to  a  certain  extent  arbitrary. 

It  cannot  be  taken  too  large  due  to  the  smallness  of  vaporization 
or  too  small,  following  the  condition  of  ordered  flow  and  hori¬ 
zontally  of  the  flow  lines.  The  photographs  of  a  vaporizing  flame 
make  it  possible  to  assume  that  0.7y*  <  X'  <  1.2y».  Let  us  write 
a  system  of  gas-dynamic  conservation  equations  for  flows  I  and  II 
between  sections  1-1  and  1»-1*.  Figure  5  depicts  the  appropriate 
monitoring  surfaces:  flow  I  -  monitoring  surface  BKSRNB ,  flow  II  - 
monitoring  surface  KAMRSK,  total  flow  I  +  II  -  monitoring  surface 
BAMNR .  The  areas  of  the  characteristic  flow  cross  sections  F',  F  , 

F*  are  known>  lf  we  know  shape  of  the  flame  and  the  shock  wave! 

It  is  possible,  for  example,  to  approximately  calculate  them  as  areas 

of  ellipses  with  semi-axes  y’/2,  z',  yQ/2,  Z(},  y */2,  z*.  The  system 

will  consist  of  equations  of  flow  rate,  energy  and  state  for  flows 

I  and  II  -  i.e.,  of  six  equations  and  of  the  momentum  equation  for 

total  flow  I  and  II  -  a  total  of  seven  equations  relative  to  eight 

unknowns  F^  W,  P»,  T»,  p»,  w",  T«,  p«  (equality  of  the  static 

pressures  P"  *  P'  over  entire  flow  cross  section  is  assumed  f  ) 

nr  * 

The  equation  of  conservation  of  flow  rate  in  flow  I  is 

^i»iPi  =  f'w'p'  ~  m,  =  m'.  ( 7 ) 

The  equation  of  the  conservation  of  energy  in  flow  I  is 

m'  ['••  +  ~2f}  +  [l*.  +  ^r]  =  [i  +  ^f\  +  ( 8 ) 

+  [»*  +  -^j-j  , 

where 

'».  =  'h  +  cPi7’1;  i«  =  iH  +  cpr  -  the  enthalpies  of  the  gas; 
im‘ =  *hm  +C*T<'  =  :'hik  +  cmT„  —  the  enthalpies  of  the  liquid  fuel; 
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i  #  i 

H*  HW 


the  initial  enthalpies.  For  the 
majority  of  real  flows  the  enthalpy 
®  liquid  is  small  in  comparison 
with  the  enthalpy  of  a  gas. 


The  gas  state  equation  in  flow  I  is 

P'lp'^gRT'. 


(9) 


The  equation  of  the  conservation  of  flow  rate  in  flow  II  +  I  is 
F4>lWl  -  m‘  +  -  F'pV  +  (F„  -  F')p'w\  (10) 


The  equation  of  conservation  of  energy  in  flow  II  is 


where  +  ejr. 


,  Auf  .  j„.i 


(ii) 


The  gas  state  equation  in  flow  II  is 

P'/p"  *•  sHT.  (12) 

The  equation  of  conservation  of  momentum  for  the  total  flow  I 
and  II  in  projections  on  the  Ox  axis  is 

Pi  Pm  +  Pi  w*Fh  +  mm  Vu  *=  P'F„  4.  mV  4- mV,  M  5 ) 

where 

m"  =  (Fu  —  F')  pv. 


This  system  of  equations  is  open;  the  number  of  equations  for 
unity  is  greater  than  the  number  of  unknowns.  Solving  it  it  is 
possible  to  express  all  unknowns  through  any  one  of  them/  As  such 

a  free  unknown,  let  us  select  the  Mach  number  M*  of  gas  in  cross 
section  F1 


M*  _  «* 

(1*0 

Assuming  for  simplicity  the  assumption  that  the  specific  heats 
°r03S  3ectl0ns  'l  «e  Identical  and  k'  -  kj,  the  system 
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off  7T2TT  (7,'<U>  Can  be  redUCed  t0  “  equation 

a  egree  relative  to  dimensionless  velocity  w'/w  ■  let 
us  call  It  an  lndiclal  equation  *  1’  let 

F& “I""-  (15) 

Its  coefficients  apart  from  the  Initial  data  and  character- 
istlc  cross  sections  of  the  flame  contain  the  free  parameter 
ecordlng  determined  value  M',  It  Is  possible  to  solve  equation 

i./T  anp,^th  ?“  °f  3l"ple  formulas  t0  find  the  other  unkno»ns 

T  /T1»  P  /P1»  w  /w\  _T"/T1#  3',  P1/P«, 

Let  us  estimate  the  unknown  parameters  -  the  velocity  v’  and 
emperature  T'  of  the  drops  In  the  cross  section,  going  Into  the 
initial  equations.  Relative  velocity  u'  should  not  exceed  the 
mum  speed  umax,  at  which  a  drop  will  break  apart  (In  cross 

c™^°t  I!  the  Pr°CeSS  °f  UqUld  «»l»*Wtlon  la  already 

completed);  v'  •  „■  .  u*.  Velocity  u  Is  determined  by  the 

fragmentation  number  [2]  or  by  the  critical  Weber  number 


Then 


because  values  p'  and  px  are  close. 


WeKp  =c 

a 

'  p'ai  ~  r  0l»,  . 


Thus,  u*  can  vary  within  the  narrow  limits  0  <  u»  <  u 

f  ii/J  j  ,  ~ ” 


magnitude  a  is  noticeably  less  than  w\  therefore  the  estimate' 
u  can  be  rather  rough.  In  the  calulatlons  let  us  take  the 
average  In  the  Indicated  range  of  values 


'a  Wo'“ 

We"p  sr  10. 


alue  Is  estimated  roughly  as  the  average  between  the 
initial  temperature  and  the  boiling  point,  since  the  enthalpy  of 
e  drops  Is  small  m  comparison  with  the  enthalpy  of  the  gas 
ar  atlon  In  u-  and  T;  has  comparatively  little  effect  on  the 
nal  results.  Thus,  In  accordance  with  everything  that  has  been 
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stated,  in  the  given  cross  section  F*  we  obtain  a  whole  set  of 
solutions  for  the  desired  values  as  a  function  of  the  free 

parameter  -  M  number.  The  question  of  selecting  any  single  solu-  ' 
tion  arises. 

A  considerable  part  of  the  flow  flowing  into  the  flame  after 

the  powerful  shock  wave  has  the  number  M  <  1.  The  supersonic  flow 

elements  (the  zones  beyond  the  section  of  the  weak  shock  wave),  in 

passing  through  dense  cluster  of  drops,  apparently,  also  become 

subsonic  except,  perhaps,  for  the  peripheral  region).  Subsequently, 

the  removal  of  work  from  the  gas  on  the  entrainment  of  drops  -  the 

process  which  is  considered  predominant  here,  accelerates  flow. 

Consequently,  the  Mach  number  of  the  subsonic  flow  on  the  average 

increases,  approaching  unity,  but  it  does  not  attain  it  (since  the 

sign  of  the  gas-dynamic  effect  does  not  change):  M  <  M  <  1. 

n  i  c 

As  a  hypothesis  let  us  suppose  that  in  the  exit  cross  section 
of  the  initial  flame  section  P'  the  M'  number  approaches  unity 
(remaining  somewhat  less  than  it),  i.e.,  M'  =  1.  Thus,  the  exit 
section  F'  is  close  to  critical.  It  can  also  be  noted  that  at 
other  comparatively  close  values  of  M'  number,  for  example,  at 

M'  *  0.75,  the  values  of  all  the  desired  parameters  do  not  vary 
very  significantly. 

The  accepted  hypothesis  of  "critical  cross  section"  at  the 
exit  from  the  initial  section  of  the  atomization  flame  closes  the 
system  of  equations  and  makes  it  possible  to  select  a  single  solu¬ 
tion.  Other  methods  of  selection  are  also  possible.1  The  table 
gives  the  results  of  the  solution,  obtained  on  a  computer,  of 
systems  of  equations  (7)-(15).  Indicial  equation  (15)  has  two  real 
very  close  roots  w'/w^  (for  each  M’).  For  purposes  of  uniformity 


ca?.be  Perfected,  by  closing  the  system  of  equa¬ 
tions  with  the  equation  of  the  gas  entropy  of  flow  II,  which  ha* 
passed  the  shock  wave  (eliminating  the  hypothesis  M'  s  l) 
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it  is  always  possible  to  take,  for  example,  the  first  root.  The 
parameters  in  the  table  show  that  the  flow  in  the  cross  section  is 
sharply  decelerated  and  is  more  heated,  than  the  incoming  flow,  and 
the  gas  flow  rate  in  the  flame  is  noticeably  less  than  the  gas  flow 
rate  in  the  incoming  flow,  which  leads  to  high  concentrations  of 
fuel  3*  and  to  small  values  of  the  coefficient  of  air  excess. 


Data  of  the  theory  and  experiment  for  gasoline 
(M  =  2.85;  P,  »  1  atm(abs . ) ;  T  ■  293°K). 
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79 
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Designations:  aTM  ■  atm(abs.);  m/cbk  ■  m/s; 
mkm  =  ym. 


A  COMPARISON  OP  THE  DATA  OP  THE 
THEORY  AND  THE  EXPERIMENT 


Having  generalized  what  has  been  stated,  we  obtain  a  final 
method  for  the  approximate  computation  of  flow  processes  with  phase 
transitions  in  the  atomization  flame. 

1.  Prom  the  initial  data  of  the  incoming  flow  and  the  fuel 
feed  (M^,  P^,  T^,  dc,  AP,  T  )  an<*  the  physical  constants  of  the 

gas  and  fuel  let  us  determine  the  flame  boundaries,  its  character¬ 
istic  cross  sections  and  the  drop  coarseness. 
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2.  Let  us  find  all  the  parameters  (w‘,  P',  T',  0')  at  inlet 
P’lnto  the  main  eectlon  of  the  flame,  having  selected,  for  example 


3.  Let  us  place  the  obtained  values  as  Initial  values  Into 
system  of  equations  (l)-(«)  and  let  us  carry  out  the  calculation 
of  the  flow  of  a  two-phase  gas-drop  mixture,  considering  the  Intra- 
flame  flow  one-dlmenslonal  and  noninteracting  with  the  surrounding 
medium.  In  this  case  It  Is  possible  to  determine  all  the  charac¬ 
teristics  -  volatility  l,  drop  velocity,  gas  parameters  -  as  a 
function  of  length  or  time. 


Let  us  carry  out  the  calculation  of  flame  volatility,  as  for 
a  flow  In  a  cylindrical  pipe,  at  the  inlet  to  which  the  two-phase 
flow  has  the  initial  parameters  w',  P',  T',  6'.  We  will  charac¬ 
terize  the  lengths  by  the  Interval  of  total  vaporization  X  obtained 
in  the  calculation,  where  all  the  liquid  particles  of  the  real 
flame  vaporize;  the  dimensions  of  the  largest  drop  a  of  the 
atomization  spectrum  are  put  Into  the  calculation.'  pressure 
in  the  isolated  flow  will  differ  somewhat  from  pressure  inside  the 
lame.  However,  In  the  range  of  obtained  values  P  this  effect  on 
e  vaporization  interval  is  comparatively  small  (8-15*) .  The 
results  of  theory  X  are  compared  with  the  data  of  experiment  L 
w  th  respect  to  injection  into  a  free  supersonic  jet  in  the  calcu¬ 
lated  gas  exhaust  heated  to  1050  and  870°K  and  number  M,  »  2.85. 

B-70  gasoline  was  fed  from  the  surface  of  a  cooled  atomizer  with 
ifferent  dfi  and  AP.  The  interval  of  total  vaporization  L  was 
determined  with  the  aid  of  photographs  obtained  by  Toepler^s 
method  (Pig.  6).  The  experiment  showed  that  the  fuel  manages  to 
vaporize  on  comparatively  short  (of  the  order  of  tens  millimeters) 
sections.  A  small  interval  of  vaporization  in  the  flow,  where  a 
drop,  it  seemed,  moves  with  supersonic  speed,  was  at  first  somewhat 
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strange.  However,  our  theory  gives  a  simple  explanation  for  It. 

The  flow  In  the  main  section  with  initial  subsonic  parameters  in 
section  F'  (their  values  w'  and  T'  are  placed  in  the  table)  is 
kept  subsonic  to  the  end  (as  a  result  of  intense  heat  withdrawal 
during  vaporization  of  the  drops).  This  leads  to  small  X  (see 
Fig.  4);  a  drop  of  the  same  dimensions  in  a  supersonic  flow  would 
nave  a  several  times  greater  vaporization  interval  (small  retention 
times ) . 

IKAPHICS  NOT  riKPHODUCri  LF 


Fig.  6.  A  vaporized  atomization  flame  of 
gasoline  in  a  supersonic  flow  ( M ^  =  2.85, 

Tj  =  1050°K,  P.  =  1  atm(abs . ) ,  AP  =  40 

atm(gage),  J0  *  1  mm). 

Taking  the  complexity  of  the  studied  process  and  the  degree 
of  approximation  of  the  solution  method  into  account  which  consist;' 
of  several  separate  stages,  the  obtained  agreement  of  the  result.-, 
of  the  theoretical  and  experimonta  1  investigations  can  be  taker;  to 
be  satisfactory. 
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